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I. INTRODUCTION 

Nonzero ground-state entropy, So ^ 0, is an important 
phenomenon in statistical mechanics. An example of this 
is water ice, for which Sn_= 0.82±0.05 cal/(K-mole), i.e., 
So/R = 0.41 ± 0.03 0 - 0 . A model with So 7 ^ 0 is the 
g-state Potts antiferromagnet (AF) on a lattice A for suf¬ 
ficiently large q @j]. This subject also has an interesting 
connection with graph theory, since the partition func¬ 
tion of the g-state Potts antiferromagnet at temperature 
T = 0 on a graph G is 

Z{G,q,T = D) PAF = P(G iq ) , (1.1) 

where P(G,q) is the chromatic polynomial of G, which is 
equal to the number of ways of coloring the vertices of G 
with q colors subject to the constraint that adjacent ver¬ 
tices must have different colors. Such a color assignment 
is called a proper vertex g-coloring of G. The minimum 
number of colors required for a proper vertex g-coloring 
of the graph G is called the chromatic number of the 
graph, denoted x(G). We will focus here on regular N- 
vertex lattice graphs Ajy and, in particular, on the ther¬ 
modynamic limit N — > 00 (with appropriate boundary 
conditions), which will be denoted simply as A. In this 
limit, the ground-state (i.e., zero-temperature) degener¬ 
acy per vertex (site) of the g-state Potts antiferromagnet 
on A is given by 

W(A,q) = lim P(A JV ,g) 1/JV , ( 1 . 2 ) 

N-¥ 00 

and the associated ground-state entropy per site is given 
by So(A, g) = ks In IF(A, g). It will be convenient to 
express our bounds on the ground-state entropy per site 
in terms of its exponent, e s °( A,q ^ kB = W(A,q). 

In @,0 by S.-H. Tsai and one of us (RS), lower bounds 
on W(A,q) were derived for the triangular ( tri ), honey¬ 
comb (he), (4-8 2 ), and sqd lattices. Here an Archimedean 
lattice A is defined as a uniform tiling of the plane 
with a set of regular polygons such that all vertices are 


*(a): Permanent address; (b) Address on sabbatical 


equivalent. Our notation for an Archimedean lattice fol¬ 
lows the standard mathematical format 0,03, namely 
A = 01 i ph ), where the product is over the regular poly¬ 
gons pi that are traversed in a circuit around a vertex 
and ca > 1 refers to possible contiguous repetitions of a 
given type of polygon in such a traversal. The sqd lat¬ 
tice is a nonplanar lattice formed from the square lattice 
by adding edges (bonds) connecting the two sets of di¬ 
agonal next-nearest-neighbor vertices in each square. In 
Q, Shrock and Tsai derived corresponding lower bounds 
on W(A,q) for all Archimedean lattices A and their pla¬ 
nar duals, using a coloring compatibility matrix (CCM) 
method employed earlier by Biggs for the square (sq) lat¬ 
tice 81, in combination with the Perron-Frobenius theo¬ 
rem [9| and a theorem giving a lower bound on the max¬ 
imal eigenvalue of a symmetric non-negative matrix [13- 

In this paper we introduce several generalizations of 
the method used in 1-0 and apply these to derive im¬ 
proved lower bounds on W(A,q) for several lattices A. 
Refs. 0 and 0 6 ] also used CCM methods to derive up¬ 
per bounds on W( A, g). However, it was shown in 0, 0 
that, while the upper bounds were moderately restrictive, 
the lower bounds were very close to the actual values of 
W(A, g). Therefore, as in [7|, we focus here on the lower 
bounds on W(A,q). 

This paper is organized as follows. In Section[n]we ex¬ 
plain the basic coloring compatibility matrix method. In 
Section Hill we discuss our generalizations of this method. 
In Sections HVlIVIl we apply our generalized methods to 
derive new and more restrictive lower bounds on W (A, g) 
for the square, triangular, and honeycomb lattices. In 
Sections IVIII and IVIIII we present corresponding results 
for two heteropolygonal Archimedean lattices, namely, 
the (4 • 8 2 ) and (3 • 6 • 3 • 6 ) (i.e., kagome) lattices. In 
Section |IX] we report results for the sqd lattice. In Sec- 
tion[X]we compare the large-g Taylor series expansions of 
our lower bounds for the various lattices with the large- 
q series expansions of the actual W functions for these 
respective lattices. Our conclusions are given in Section 
m We list some results on r-partite lattices in Appendix 
S the lower bounds on W (A, q) for Archimedean lattices 
A from 1-0 in Appendix [Bl and some higher-degree al¬ 
gebraic equations that are used in the text in Appendix 
0 
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II. BASIC CALCULATIONAL METHOD 

In this section we explain the basic calculational 
method used in 1-1 to derive lower bounds on W (A, q). 
In the next section we generalize this method in several 
ways. We consider a sequence of (regular) lattices of type 
A of length L x = n vertices in the longitudinal direction 
and width L y = m vertices in the transverse direction. 
In the thermodynamic limit n —>■ oo, m —> oo with the 
aspect ratio m/n finite, the boundary conditions do not 
affect W(A,q). It will be convenient to take periodic 
boundary conditions (PBCs) in both directions. If a lat¬ 
tice A is r-partite, then to and n are chosen so as to 
maintain this property. 

The construction of the coloring compatibility matrix 
T begins by considering an n- vertex path V n in the lon¬ 
gitudinal direction on A. The number of proper vertex 
q-colorings of V n is the chromatic polynomial P(V n ,q). 
Now focus on two adjacent parallel paths, V n and V' n . 
Define compatible proper q-colorings of the vertices of 
these adjacent paths as proper g-colorings such that no 
two adjacent vertices on V n and V' n have the same color. 
One can then associate with this pair of adjacent paths 
an M x Af dimensional symmetric matrix T, where Af = 
P(Pn,q) = P(n,9), entries Tp„,-p; = T-p;,-p„ = 1 
or 0 if the proper g-colorings of V n and V' n are or are not 
compatible, respectively. This matrix is thus defined in 
the space of allowed color configurations for these adja¬ 
cent paths. 

It follows that, for fixed m and ?r, 

P(A mX n, <?) = Tr(T m ) . (2.1) 

For a given n, since T is a nonnegative matrix, one can 
apply the Perron-Frobenius theorem [9] to conclude that 
T has a real positive maximal eigenvalue A max . Hence, 
for fixed n, 

lim Tr(T m )™" = (A max )" . (2.2) 

m—too 

Therefore, taking the n —> oo limit, 

W(A,q) = lim (A max )™ . (2.3) 

n—>• oo 

Let us denote the column sum 

Af 

K j (T)=^T ij , (2.4) 

i=1 

which is equal to the row sum, 


Note that S(T)/Af is the average row sum (equal to the 
average column sum). 

For a general nonnegative Af x Af matrix A, Q, one 
has 

min{Kj(A)} < A max (A) < max{ft,,(A)} (2-7) 

and 

min {pj{A)} < A max (A) < max{pj(A)} (2.8) 

for j = Since T t = T, these are equivalent here. 

One also has the following more restrictive one-parameter 
family of lower bounds depending on the parameter k , for 
a symmetric nonnegative matrix T 

- JJ 1 <A maxiT). (2.9) 

Refs. 1-1 derived lower and upper bounds on 
W(A,q) using the k = 1 special case of (12.91) . We will 
denote a generic lower bound on W (A, q) with the sub¬ 
script £, as W(h.,q)i. We will distinguish specific lower 
bounds that we obtain with the additional subscripts b 
and k , as explained below. The lower bounds obtained 
in 1-1 were for 6=1 and k = 1. Refs. 1-0 stud¬ 
ied how close the upper and lower bounds obtained on 
W (A, q) were to the actual values of W (A, q) for a num¬ 
ber of lattices, where the latter were determined mainly 
from Monte Carlo calculations, augmented by large-g se¬ 
ries expansions together with a few exact results. It was 
found that for a given lattice A, as q increases beyond 
the region of x(A), the lower bounds rapidly approach 
very close to the actual value of W (A, q). 

We next introduce some notation that will be used be¬ 
low for reduced functions obtained from W (A, q) which 
will be analyzed in the large- q limit. This large- q limit 
is the natural one to consider for chromatic polynomials, 
since the constraint in a proper g-coloring of the vertices 
of a graph, namely that no two adjacent vertices have 
the same color, becomes progressively less restrictive as 
the number of colors increases to large values. The chro¬ 
matic polynomial of an arbitrary A-vertex graph G is a 
polynomial of degree N, and consequently, W (A, q) ~ q 
as q —> oo. In order to deal with a finite quantity in 
the q —> oo limit, one therefore considers the reduced (r) 
function W r (A,q) = W(A,q)/q. A variable equivalent to 
1 /q that is convenient to use for a large-g series expansion 
of W r (A,q) is y = l/(g— 1). These large-g (i.e., small-y) 
series expansions are normally given for the function 


N 

Pj (T) = J2Tji, ( 2 . 5 ) 

i— 1 

(since T t = T) and the sum of all entries of T as 

M 

Sir) = Tij . (2.6) 

i,j =1 


W( A,y) 


W r (A,q) 

(1 - q~ !) a a/2 


W(A,q) 

qil-q- 1 )^/' 2 


( 2 . 10 ) 


where Aa is the lattice coordination number of the lat¬ 
tice A (i.e., the degree of the vertices of A). In terms of 
the expansion variable j/, these series thus have the form 
W(A,y) = 1 + Analogously, for the ex¬ 

pansion of our lower bound, we define the reduced lower 
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bound function W(A,y)t as 


W(A,y) t 


W(A,g) t 

q( 1 — g _1 ) A a/2 


( 2 . 11 ) 


Before proceeding, we note a subtlety in the definition 
of W(A,q). As pointed out in [12], the formal eq. (11.21) 
is not, in general, adequate to define W( A, q) because of 
a noncomnrutativity of limits 

lim lim P{A N ,q)^ N ^ lim lim P(A N ,q) 1/N 

iV—>-oo <?—>-g s q—*q s N—>oo 

( 2 . 12 ) 

at certain special points q s . We denote the definitions 
based on the first and second orders of limits in (12.121) as 
W (A, q)o Nq and W (A, q)jj qN , respectively. This noncom¬ 
mutativity can occur for q < q c { A), where q c { A) denotes 
the maximal (finite) real value of q where W (A, q) is non- 
analytic [12]. These values include q c (sq ) = 3, q c (tri) = 
4, and the formal value q c {hc) = (3 + \/5)/2 = 2.618... 
[j, [HI for the square, triangular, and honeycomb lat¬ 
tices. As explained in [12], the underlying reason for 
the nonconrmutativity is that as q decreases from large 
values, there is a change in the analytic expression for 
W(A,q) as q decreases through the value q c ( A). We do 
not have to deal with this complication here because ele¬ 
mentary results yield exact values of W(sq, 2), W{hc , 2), 
W((4- 8 2 ), 2), and W(tri, 3) (see Eqs. (1AH) and (1A2ID . 
namely 

W{sq, 2) = W(hc, 2) = W((4-8 2 ), 2) = 1, W{tri, 3) = 1. 

(2.13) 

Hence, our lower bounds are not needed at the respective 
values q = 2 for the square, honeycomb, and (4 • 8 2 ) 
lattices or for q = 3 on the triangular lattice, and we 
therefore focus on their application to q > 3 for A = 
sq, he, (4 • 8 2 ) and to q > 4 for A = tri, and similarly 
for other lattices. 


III. GENERALIZED COLORING 
COMPATIBILITY MATRIX METHOD 

A. Coloring Compatibility Matrix Joining 
Adjacent Strips of Width b 

The lower bounds on W( A, q) derived in |-1 for vari¬ 
ous lattices A used Eq. m with T being a coloring com¬ 
patibility matrix joining adjacent paths and with k = 1. 
Here we generalize this method in several ways. Our first 
generalization is to use a coloring compatibility matrix 
that joins adjacent strips of width b > 2 vertices, rather 
than adjacent one-dimensional (b = 1) paths. For sim¬ 
plicity, we explain this for the square lattice; similar dis¬ 
cussions apply for other lattices. We define the matrix T 
to enumerate compatible colorings of a strip of transverse 
width b vertices and an adjacent parallel strip of width 
b and arbitrary length L x vertices, with cyclic boundary 
conditions. (Here, by cyclic boundary conditions for a 


given strip, we mean in the x, i.e., longitudinal, direction 
along this strip). The condition that these strips are ad¬ 
jacent is equivalent to the statement that they share a 
common set of edges. Thus, this CCM is an Af x Af ma¬ 
trix, where M is the chromatic polynomial for the cyclic 
strip of width b vertices and arbitrary length L x , with 
cyclic boundary conditions. For this CCM, the sum of 
elements S(T) is equal to the chromatic polynomial of a 
strip of width L y = 2b — 1 vertices and arbitrary length 
L x vertices with cyclic boundary conditions. These chro¬ 
matic polynomials of lattice strips of a fixed width L y 
and arbitrarily great length L x with periodic boundary 
conditions in the longitudinal direction and free bound¬ 
ary conditions in the transverse direction have the form 

Ly np(Ly ,d) 

P( A, Ly x L x , cycl, g) = ^ c(d) ( x ^,L v ,dj) Lx 

d =0 j=1 

(3.1) 

with 

c (d) =£(-l y ( 2d ~ J )q d - j , (3.2) 

3=0 \ 3 J 

where (£) = a\/[b\{a — 6)!] is the binomial coefficient. 
For a table of the np(L y ,d), see [Hj]. Because of the 
limits m and m, only the largest Aa ,L v ,d,j enters in 
the lower bound (12.91) in the thermodynamic limit. As 
specific studies such as 0-E3 showed, the dominant A 
for the values of q of relevance here is Aa,l„,o,i- 

Applying this generalization of the coloring compati¬ 
bility matrix in combination with the k = 1 case of m , 
we derive the new lower bound for b > 2: 


W(A,q)>W(A,q) e , b ,i , 

where 

Aa,26—1,0,1 
Aa, 6 ,o,i 


W(A,q)f. bt i = 


(3.3) 


(3.4) 


The final subscript, 1, in W(A, q)e-b,i in (13.31) and (13.41) 
is the value of k. 

The corresponding lower bound for W (A, y) is 

W(A,y)>W(A,y) eibtl , (3.5) 

where, in accordance with Eq. (12.101) . 


W(A,y) e , btl = 


W(A,q) e , btl 

g(l-g-!) A A/2 ’ 


(3.6) 


with Aa being the coordination number of the lattice 
A, as before. The inequality (13.31) with (13.41) is actually 
an infinite family of lower bounds depending on the strip 
width b = 1,2,..., and similarly with (13.51) and (13.61) . This 
is one of our two major results, which we will proceed to 
apply to a number of different lattices. The special case 
5=1 was previously used in [g] and 1-0 to derive 
lower bounds which we denote here as W{A, g)^ ; i,i and 
correspondingly W(A, y)^ ; i,i. Our generalization in this 
subsection is to 5 > 2 with k = 1. 
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B. Coloring Compatibility Matrix Acting k Times 
Joining Paths of Width b = 1 


Our second generalization is to use a coloring com¬ 
patibility matrix method that involves paths (i.e., one¬ 
dimensional strips, with b = 1) on A that are separated 
by k edges, where k > 2 , rather than the situation with 
b = 1 and k = 1 considered in 1 - 0 . where the paths 
were adjacent. This means using the coloring compat¬ 
ibility matrix T defined as connecting adjacent paths, 
and having it operate k times, with k > 2. Hence, 
M = P(C Lx ,q) and S(T k ) is the chromatic polynomial of 
a strip of width L y = k + 1 vertices and arbitrary length 
L x vertices with cyclic boundary conditions. Again, only 
the dominant Aa ,L y ,d,j terms enter in (12.91) in the thermo¬ 
dynamic limit. Using this method in combination with 
m . we derive the lower bound 


W(A,q) > W(A, q)e-,l,k , 


where 




AA,fc+l,0,l 
Aa, i,o,i 


1 

T 


(3.7) 


(3.8) 


In (13.71) and (13.81) . the first subscript after l\ is b = 1. 

An important theorem extending the result (EH) is 
that for a symmetric nonnegative matrix T fl 8 | . 



AT 


is an increasing function of k . 


(3.9) 


It follows that, for the physical range of q of relevance for 
our application to a lattice A, 


C. Measures of Improvement of Bounds 

For a lattice A and a given q , we define the ratio of a 
lower bound W (A, q)e.-b,k to the actual value of W{A,q) 
as 


Pa ,q]£;b,k 


W (A, q)t-b,k 
W(A,q) 


(3.13) 


This ratio is useful as a measure of how close a particular 
lower bound W(A,q)e-b,k is to the actual value of the 
ground-state degeneracy per vertex, W(A,q). For most 
lattices and values of q, the value of W (A, q) is not known 
exactly, but rather is determined for moderate values of 
q by Monte Carlo simulations, as discussed in [t| [f$ and, 
for larger values of q, by large-g series expansions [ 20 j . 
Special cases of A and q for which exact results are known 
will be noted below. 

An important property of our new lower bounds is 
that, for a given lattice A, they are larger than and hence 
more restrictive than the bounds VF(A, q)e-i i derived in 
1 - 1 - Since the lower bounds W(A, q)e-,i,i were very 
close to the actual values of W (A, q) for all but the lowest 
values of g, our improved lower bounds are even closer to 
these actual values. For the same reason, our new lower 
bounds yield the greatest fractional improvement for low 
to moderate values of q and are only slightly greater than 
W (A, q)e-i,i for larger values of q. This will be evident in 
our explicit results. For our present discussion, we take 
T to be the matrix that acts k times mapping a strip 
of width b to an adjacent strip of width b on A. Then 
the theorem (13.91) and its corollary (13.101) imply that, for 
fixed b , the ratio of our lower bound W(A, q)e-b,k to the 
actual value W(A,q) is an increasing function of fc, i.e., 


RA,q-/;b,k is an increasing function of k . (3.14) 


W(A,q)e- t i,k is an increasing function of k . (3.10) 

The corresponding lower bound for W (A, y) is 

W{A,y) >W(A,y) t;lik , (3.11) 

where, in accordance with Eq. (12.101) . 


W(A,y) e . xk = 


W(A,q)e-i t k 

q(l — q~ 1 ) A ^/ 2 ' 


(3.12) 


Again, the inequality (13.71) with (13.81) is actually a one- 
parameter family of lower bounds depending on the pa¬ 
rameter k = 1,2,..., and similarly with (13.111) and (13.121) . 
This is the second of our major results. The special case 
k = 1 (with b = 1 ) was previously used in the 

generalization presented in this subsection is to k > 2 
with 6=1. We have also carried out further generaliza¬ 
tions of lower bounds on W (A, q) with both b > 2 and 
k > 2. These are more complicated and will be presented 
elsewhere. 


That is, as k increases, the lower bound W (A, q)t,b,k be¬ 
comes more restrictive. From our analysis, we also find 
that for fixed k = 1 and b > 2 , 

RA,q-/\b,i is an increasing function of b . (3.15) 

For a given A and q, it is also of interest to compare the 
various lower bounds with each other. For this purpose, 
we define the ratio 

D _ W(A,q) ( . b ,k 

R A, q ;(b,k)nb',k')~ W{Kq)t . h , k , ' ( 3 ‘ 16 ) 

By the same argument, theorem (13.91) and its corollary 
(IB. 101) imply that for a given lattice A, our new lower 
bounds W(A,q)t-\^k improve on the bound W(A : q)i-i t i 
derived in |l]-[ 8 ]: W(A,q) e .,i,k > TU(A, < 7 )^ 1 ,i, i.e., 

7?A, g; (i,fe)/(i,i) >1 for k > 2 . (3.17) 

We observe also that 

#A,g;(&,i)/(i,i) >1 for b > 2 . (3.18) 
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As will be evident from our explicit results, for the range 
of q that we consider, these inequalities are realized as 
strict inequalities. As noted above, since the latter lower 
bounds W(A, q)e-,i,i are very close to the actual values 
of W(A,q), even for q only moderately above y(A), as 
shown in Table I of [H] and Tables I-III of [f|, our new 
bounds are even closer to these actual values of W( A, q). 
In all cases, we find that the ratios approach unity rapidly 
in the limit q —> oo. 

A major result of Ref. [?] was the derivation of 
general formulas for the lower bound W(A,q)^.ii and 
W(A, y)qi,i for all Archimedean lattices and their (pla¬ 
nar) duals (Eqs. (4.11), (4.13), (5.1), and (5.2) in Q). As 
will be evident below, aside from the basic theorems, our 
new lower bounds W(A,q)t,b,k with 5 > 2 and/or k > 2 
do not have such simple general formulas. However, as 
noted, they do provide a useful improvement on the ear¬ 
lier W(A, q)e-i,i lower bounds, especially for q values not 
too much larger than x(A). 


In terms of W(sq,y), given by (12.101) with A = sq and 
A = 4, the lower bound is the 5=1 case of (ESI) with 
A = sq, namely 

W(sq,y)t ]1A = l + y 3 , (4.3) 

as listed in Table III of Q. 

A. CCM Method with b = 2, 3 and k = 1 

We first use our generalized method with the color¬ 
ing compatibility matrix relating the allowed colorings of 
a width 5 = 2 cyclic ladder strip of the square lattice 
to those of the adjacent 5 = 2 strip. For this, we need 
the dominant term in the chromatic polynomial for the 
square-lattice strip of width 25 — 1 = 3 for the relevant 
range of q > 3. This chromatic polynomial was calcu¬ 
lated in [l4|, and the dominant term is Agg^o,!, namely 


IV. SQUARE LATTICE 


A 


_ 1 

sg,3,0,l — ^ 


(< Z - 2)( 9 2 


3 g + 5) 


As noted above, since the value W(sq, 2) = 1 is known 
exactly by elementary methods, we focus on the applica¬ 
tion of our new lower bounds to the range q > 3. We first 
recall the result for the case 5 = 1, k = 1. With T being 
the coloring matrix connecting adjacent rows or columns 
of a square lattice, and with the application of the k = 1 
special case of the theorem (12.91) . one has 

W(sq,q) > W(sq,q) e -i,i , (4.1) 

where [§| 

W(sq,q)e-i,i = - - — + ’* ■ (4.2) 

q- 1 

I 


(q 2 -5 q + 7)(q 4 - 5 q 3 + 11 q 2 - 12q + 8) 


1/2 


(4.4) 


This term is also the dominant A in the chromatic poly¬ 
nomial for the strip of the square lattice with transverse 
width L y = 3 vertices and arbitary length, with free lon¬ 
gitudinal and transverse boundary conditions [T9] . Our 
lower bound with 5 = 2 (and k = 1) then reads 

W(sq, q) >W(sq, q)e - 2 ,i , (4.5) 

where (with A sg , 2 , 0,1 = q 2 — 3g + 3) [2l| 


W(sq,q)e- 2 ,i = 


A S( j,3,0,1 


(q - 2 )(q 2 - 3q + 5) + (q 2 -5 q + 7)(q 4 - 5 q 3 + 11 q 2 - 12 q + 8) 


1/2 


A 


sq, 2,0,1 


2 (q 2 - 3q + 3) 


(4.6) 


Using the analytic results (ira and (ira . we have 
proved the following inequality (for q > 3): 

W(sq,q)(- 2 ,i>W(sq,q)e-i,i (4.7) 

In terms of the ratio R S q.q;( 2 , 1 )/(i,i) ? 

-Rsg,< 7 ;( 2 ,l)/( 1 , 1 ) ^ 1 (4.8) 

The inequality (14.81) means that our new lower bound, 
63, is more stringent than the previous lower bound 
63 obtained with the CCM method with 5 = 1 and 
k= 1. 


Ref. Q showed that as q increases, R sq , q - 1,1 rapidly 
approaches extremely close to unity. For example, for 
q = 4, 5, 6, R sq , q -i,i is equal to 0.9984, 0.9997, and 
0.9999 (see Table I in 1), respectively, and it increases 
monotonically with larger q. Our improved lower bound 
63 on W(sq, q) is therefore even closer to the respective 
actual values of W(sq, q). As will be discussed below, this 
is also true of our other new lower bounds using 5=1 
and k > 2. We note that if one were formally to extend 
the range of applicability of (14.71) down to q = 2, it would 
be realized as an equality, and if one were to extend the 
range of applicability of (14.81) to 2 < q < 00 , it would be 
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realized as an equality at q = 2 and in the limit q —> oo. with (14.61) improves on this. For q = 3, we have 
For the previous lower bound W(sq , <?)qi,i, the largest 
deviation from the actual value occurs at q = 3. It hap- g _l_ ^jyj 

pens that for q = 3, W(sq, 3) is known exactly [Hf: T,r '^ 


W(sq 1 3) = —== = 1.5396007... 
3 d / 2 


For the old bound, 

W(sq, 

so that 


_ 3 

q— 3 2 


W(sg, 3)^1,] 


g=3 


fF(s g ,3) 


3 5 / 2 

16 


= 0.974279 


: W(sq, 3)^2,i 

(4.9) g0 ^hat 


(4.10) 


(4.11) 


9=3 


6 


= 1.5205176..., (4.12) 


W(sq, 3)*;2,] 


i? 


S9,3;^;2,l = 


? =3 V3(5 + vT7) 


W(sg,3) 


16 


= 0.987605... 
(4.13) 


We show these ratios R sq , 3 -e-i,i and R sq , 3 - 1 - 2 ,i in Table UJ 
In terms of the function W ( sq , y ), our lower bound 
(S3 reads 


W{sq, y) > W(sq,y)t, 2 , 


l ) 


to the indicated floating point accuracy. As guaranteed 

by the general inequality (14.81) . our lower bound 63 where 


(4.14) 


(i + y) 


W{sq,y)t- 2 ,i = 


(1 - y)( 1 - y + 3y 2 ) + (1-3 y + 3y 2 )(l - y + 2y 2 - y 3 + 3y 4 ) 


1/2 


2(1 -y + y 2 ) 

We have also calculated the lower bound W{sq, q)e-b,i for b = 3, and we list the ratio R sq , 3 -/-, 3 ,i in Table |H 


(4.15) 


B. CCM Method with 6=1 and 2 < k < 5 


Next, we apply our second generalized method to the square lattice. For k = 2, our lower bound obtained using 
this method is (13.71) with (13.81) . namely 


W(sq,q) > W{sq,q)e- lt2 , 


(4.16) 


where 


W(sq,q)e- 1 , 2 = 

The corresponding lower bound on W{sq,y) is 

W{sq,y ) > W(sq,y)t. 1<2 , 

where 


A sg ,3,0,1 

1/2 

" (<Z — 2) (y 2 — 3y + 5) + 

(y 2 -5 q + 7)(y 4 - 5 q 3 + 11 y 2 - 12y + 8) 


- ^C,0,1 . 



2(9-1) 



1/2 


W(sq,y)i. 1,2 = -j=(\ + y) 


(1 - y){ 1 - V + 3y 2 ) + (1 - 3y + 3y 2 )(l - y + 2y 2 - y 3 + 3y 4 ) 

I- 


1/2 


1/2 


. (4.17) 


(4.18) 


(4.19) 


For 6=1 and k = 3, we need the dominant A in 
the chromatic polynomial for the cyclic square-lattice 
strip of width L y = 4 vertices and arbitrary length L x , 
namely A^^o,!; which was calculated in [23j (and is the 


same as the dominant A in the chromatic polynomial 
of the free square-lattice strip of width L y = 4 El)- 
This term A^^ou is the largest (real) root of the cu¬ 
bic equation (IC1I) in Appendix [Cj Our bound is then 
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W(sq,q) > W(sq,q)e- 1 , 3 , where 


W(sq,q) e - lt 3 = 


Asg,4,0,l 

9-1 



(4.20) 


In a similar manner, for 6 = 1 and k = 4, we have ob¬ 
tained the bound W{sq,q) > W(sq,q)e- 1 , 4 , where 


W(sg,g)^ ; i,4 = 


Asg,5,0,l 

9-1 



(4.21) 


and is the largest (real) root of an algebraic equa¬ 

tion of degree 7. 

As a special case of our general result (13.101) , we have 


Rsq,q;(l,4)/(1,1) > -Rsqr,qr;(1,3)/(1,1) > Rsq,q;(l,2)/(1,1) > 1 • 

(4.22) 

In the range q > 3 under consideration here, we find that 
each > is realized as >, i.e., a strict inequality. 

It is also of interest to compare our various lower 
bounds W(sq,q)e-b,i and W(sq, q)e-,i,k with each other. 
For the first two above the old case 6 = 1, k = 1, we find 


^sg,g;(2,l)/(l,2) > 1 ■ (4.23) 

That is, our lower bound with (6, k) = (2,1) is larger, 
and hence more restrictive, than our lower bound with 
(6, k) = (1,2). In the limit q —> oo, the ratio (14.231) 
approaches 1. 


C. Plots 

In Fig. [l]we plot the ratios R sq , q -(b,i)/(i,i) for 6=2 and 
6 = 3 as functions of q in the range 3 < q < 6, and in Fig. 
[2]we plot the ratios R sq , q -,(i,k)/(i,i) f° r & = 2 up to k = 5, 
as functions of q in the same range. (Here and below, 
such plots entail a continuation of the relevant expres¬ 
sions from integral q to real q.) These plots illustrate the 
result that we have proved in general, that, for a given q, 
RA,q;(i,k)/(i,i) is an increasing function of k, and also our 
result that R sq , q -,( 3 ,i)/(i,i) > R sq , q -( 2 ,i)/(i,i)- (If formally 
continued below q = 3 to q = 2, the curves reach maxima 
and then decrease; for example, R SI j, g; ( 2 ,i)/(i,i) reaches a 
maximum of 1.06 at q ~ 2.29 and then decreases to 1 as 
q \ 2, while R sq , q -(i, 2 )/(i,i) reaches a maximum of 1.03 
at q ~ 2.29 and then decreases to 1 as q \ 2.) As the re¬ 
sults in these figures show, our new lower bounds improve 
most on the earlier W ( sq , q)e-,i,i in the region of q ^ 3; as 
q increases beyond this region, the new bounds approach 
the earlier one. This feature will be evident from the 
large-g (small-?/) expansions, since the new bound and 
the earlier one coincide in the terms of the small-?/ ex¬ 
pansion up to 0(y 6 ). We also find this type of behavior 
for the new lower bounds that we have derived for other 
lattices; that is, the degree of improvement is greatest 
for the region of moderate q slightly above x(A). On a 
given lattice A, for larger q , our new bounds rapidly ap¬ 
proach the earlier one with k = 1 and 6=1; i.e., the 
ratio RA, q ;(b,k)/(i,i) rapidly approaches unity. 


1.021 
1.018: 

1.016: 

1.014) 

1.012 

- 

: 1 \ 

ratios 1 . 0 T 

1.008 \ 

: \ 

\ 

1.006: 

1.004) 

1 . 002 ) 

1 3 3.5 4 4.5 5 5.5 6 
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FIG. 1: Plot of the ratios R sq ,q-,(2,i)/(i,i) (lower curve) and 

Rsq,q-,( 3 , 1 )/(l,l) (upper curve) as functions of q for 3 < q < 6 . 



FIG. 2: Plot of the ratios R sq ,q-,(i,k)/(i,i ) f° r & = 2 to k = 5 as 
functions of q for 3 < q < 6 . From bottom to top, the curves refer 
to k = 2, k = 3, k = 4, and k = 5, respectively. 


Combining these results with the results in Table I in 
[5] and Table I in (g|, it follows that as q increases above 
the interval of q = 3 and q = 4, these lower bounds 
approach extremely close to the actual respective values 
of W(sq, q). As was evident from these tables in 11 , in 
the range q > 3, the greatest deviation of the lower bound 
W(sq,q)e-i^i from the actual value of W(sq,q) occurs at 























q = 3. It is thus of interest to determine how much 
closer our improved lower bounds are to W(sq , 3). From 
our general expression for W(sq , q)t, 1 , 2 , we calculate the 
q = 3 value 


W(sq, 3)£;i, 2 


9—3 


2 


1.510223959... (4.24) 


so that 


lF(sg,3)£;i, 2 


1 ?. 


sg,3;^;l,2 = 


g=3 


W(sq, 3) 


3 ^ 3(5 + 717 ) 


16 


= 0.9809192... (4.25) 


This ratio and the other ones discussed here are listed in 
Table |H 


V. TRIANGULAR LATTICE 

A. b = 1, k = 1 


expression [26]. For example, for q = 5, 6, 7, Rtri,q; 1,1 
is equal to 0.9938, 0.9988, and 0.9996, respectively, and 
it increases monotonically with larger q. Since our new 
lower bounds on W ( tri , q) are more restrictive than (15.11) . 
they are therefore even closer to the respective actual val¬ 
ues of W(tri, q). 

The corresponding lower bound on W(tri,y) is 
W(tri,y) > W(tri,y)e-i t i, where 

W(tri,y) e . 1A = (l-y 2 ) 2 (5.2) 

(see Table III in S). 


B. b = 2, 3, k= 1 

Here we derive a new lower bound on W(tri , q) using 
our first generalization of the CCM method with 6 = 2, 
k = 1. For this purpose, we need the chromatic polyno¬ 
mial of the cyclic strip of the triangular lattice of width 
L y = 3 vertices and arbitrary length, L x . This was cal¬ 
culated in [25|. The dominant A in (13.11) is 


Since W(tri, 3) = 1 is exactly known, we will restrict 
our consideration of lower bounds to the range q > 4. 
We recall that for 6 = 1 and k = 1, one has the lower 
bound i-0 W(tri,q) > W(tri,q)(-i t i, where 


3,0,1 — 0 


q 3 - 7 q 2 + 18 q - 17 


q b - 14<f + 81g 4 - 250 q 6 + 442 q 2 - 436 q + 193 


1/2 


W(tri,q)e-i,i = ^ • (5.1) 

As was discussed in Q, q increases beyond the lowest 
values above y(tri) = 3, this lower bound rapidly ap¬ 
proaches the known value of W(tri,q) (see Table I in 
[5| ), where the latter was determined by a numerical eval¬ 
uation of an integral representation and infinite product 


(5„ 

Combining this with A tr j, 2 ,o,i = {q — 2) 2 , we derive the 
lower bound 


W(tri,q) > W(tri,q)e- 2 ,i , (5.4) 


where 


W(tri,q) e . 2 ,i = 


^tri, 3,0,1 


q 3 - 7q 2 + 18 q - 17 + 


q 6 - 14 q 5 + 81 q A - 250 q 3 + 442 q 2 - 436g + 193 


1/2 


1 • (5-5) 


A^r-i,2,0,l 2((^ 2) 2 

The reduced function W(tri,y) is given by Eq. (12.101) with A = tri and A = 6. The corresponding lower bound is 

W(tri, y) > W(tri,y)i- 2 ,1 , (5.6) 


where 


(i + y ) 2 


W(tri , y)t, 2,1 = 


1 — 4 y + 7 y 2 — 5 y 3 + 


1-8 y + 26 y 2 - 46 y 3 + 53y 4 - 42 y 5 + 17y 6 


1 1/2 


2(1 -y? 

1 - 


(5.7) 


Our new lower bound W(tri, q)e-, 2 , 1 is larger than, and W(tri, q)t,\,i- That is, from the analytic forms (15.21) and 
hence more restrictive than the previous lower bound, 
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(EH), we have proved that (for q > 4) 


so 


Rtri,q-,{ 2,1)/(1,1) > 1 • (5-8) 

This ratio approaches 1 as q —> oo. 

As was evident in Table I in Q, the deviation of 
W(tri, q)t-i,i from the actual value of W(tri,q) was 
greatest for q = 4. Hence, it is of interest to determine 
how much closer our new lower bound W ( tri, q)e- 2 ,i is to 
the W(tri, q) for this value, q = 4. A closed-form inte¬ 
gral representation has been given for W(tri, q) [2(|; in 
particular, an explicit result is the value for <7 = 4: 


R, 


•triA:i:2,l — 


(7 + yi7)v / 3r(2/3) 3 


167T 


= 0.951669845... 

(5.13) 

We have also calculated the lower bound W(tri, q)t,b ,l 
for 6 = 3 and evaluated this for <7 = 4. For reference, we 
list the various ratios Rtri,i-,t,b,k in Table ED We see that 
W(tri, 4) 2,1 and W(tri, 4)^3^ are closer to the exact 
value of W(tri, 4) than W(tri, 4)i-,i,i- 


N 3T(l/3) 3 2tt 

W(tri, 4 = —W- = -=-= 1.460998486... 

V ' 4tt 2 v / 3T(2/3) 3 

(5-9) 

where the equivalence follows from the relation 
r(,z)r(l — z) = 7r/ sin(7rz) for the Euler Gamma function. 
We recall that 


C. b = 1, 2 < k < 5 
By the same means as above, we derive 
W(tri, q) > W(tri,q)i-i ,2 , 


(5.14) 


W(tri, 4) /;lil =| (5.10) 

SO 

2rf2/3') 3 

Rtriwxi = —r=-^- = 0.9126178746... (5.11) 

V3 7T 

(see Table I of 0)- The value of our new lower bound at 
q = 4 is 


with 


W(tri,q) e - 1,2 = 


A, 


tri, 3,0,1 


1/2 


q-1 


(5.15) 


where Atri, 3 ,o,i was given in Eq. EH- Equivalently, 

W(tri,y ) > W(tri, y)e, 1,2 , (5.16) 


W (tri, 4) e . 2 ,i = 


7 + Vri 

8 


1.3903882... 


(5.12) where 


W(tri, y)e -,!,2 = -j= (1 + y) 2 


1 — 4j/ + 7y 2 — 5y 3 + 


1-8 y + 26 y 2 - 46y 3 + 53 y A - 42 y 5 + 17 y 6 

I- 


- 11/2 


1/2 


(5.17) 


For 6 = 1, k = 3, we need the dominant A in the 
chromatic polynomial for the cyclic strip of the triangu¬ 
lar lattice of width L y = k + 1 = 4, namely, Xtri, 4 .o,i- 
This chromatic polynomial was calculated in [25| . and 
the dominant A is given as the largest root of the quartic 
equation (IC2I) in Appendix[Cj This is also the dominant A 
in the chromatic polynomial of the free strip of the trian¬ 
gular lattice with width L y = 4 and arbitrary length 0- 
We have also calculated W(tri, q)t,i,k for k = 4, 5. For 
reference, we list the various ratios Rtri, 4 -,e-,i,k i n Table ITT1 


D. Plots 


square lattice, these plots illustrate the result that we 
have proved in general, that, for a given q, J?A,q;(i,fc)/(i,i) 
is an increasing function of k, and also our result that 
Rtri,q;( 3,1)/(1,1) ^ Rtri,q-,( 2,1)/(1,1) - 


VI. HONEYCOMB LATTICE 

Since W(hc, 2) = 1 is exactly known, we restrict our 
consideration of lower bounds for the honeycomb lattice 
to the range q > 3. We recall that for 6 = 1 and k = 1, 
one has the lower bound W(hc, q) > W(hc, q)e-i i, where 

a 


In Fig. 0 we plot the ratios Rtri,q;(b, i)/(i,i) f° r 6 = 2 
and 6 = 3 as functions of q in the range 4 < q < 6, and 
in Fig. |3]we plot the ratios Rtri,q;(i,k)/(i,i) for k = 2 up 
to k = 5, as functions of q in same range. As with the 


W(hc,q) e -i,i = 


(AQ 1/2 

<7-1 


(i q 4 - 5<7 3 + 10<7 2 - 10<7 + 5) 1 / 2 


<7-1 


( 6 . 1 ) 
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FIG. 3: Plot of the ratios Rtri,q,( 2,l)/(l,i) (lower curve) and 

l^tri,q;(3,l)/(i,i) (upper curve) as functions of q for 4 < q < 6 . 



equal to 0.99898, 0.99985, and 0.99996, respectively, and 
it increases monotonically with larger q. Since our new 
lower bounds on W(he, q) are more restrictive than (16. 111 , 
they are therefore even closer to the respective actual 
values of W(hc,q). 

The corresponding lower bound on W (he, y) is 
W(hc, y) > W(hc,y)e- 1A , where 0 

W(hc,y) e . X i = (l+y 5 ) 1/2 ( 6 . 2 ) 

(see Table III in 0 )- 

For the calculation of W(hc,q)e- 2 , 1 , we need the chro¬ 
matic polynomial of the cyclic strip of the honeycomb 
lattice of width L y = 26 — 1 = 3 vertices and arbitrary 
length, L x , in particular, the dominant A. This Xhc, 3 ,o,i 
is the largest (real) root of the cubic equation (1C3I) in 
Appendix [O 2j|. This dominant A is also the input that 
we need for the calculation of W(hc, q)e-,i, 2 , since the lat¬ 
ter requires the same chromatic polynomial of the cyclic 
strip of the honeycomb lattice of width L y = k + 1 = 3 
vertices and arbitrary length, L x , in particular, the dom¬ 
inant term. This A is also the dominant term in the 
chromatic polynomial of the strip of the honeycomb lat¬ 
tice of width L y = 3 vertices and arbitrary length, with 
free boundary conditions 0 - 


VII. 4 8 2 LATTICE 


Using the CCM method with 6=1 and k = 1, Ref. 0 
derived the lower bound W( (4- 8 2 ), q) > lU((4-8 2 ), < 7 )^ 14 , 
where 


W((4-8 2 )),q) e . hl = 


(£ 4 D 8 ) 1/4 

<7-1 


(7.1) 


Equivalently, W((4 • 8 2 ),y) > W((4 • 8 2 ),y)qi,i, where 
W ((4 • 8 2 ), y) t , 1,1 = [(1 + y 3 )(l + y 7 )] 1,A . (7.2) 


We have obtained the slightly more restrictive lower 
bound IU((4 • 8 2 ),q) > W((4 • 8 2 ), <?)qi, 2 , where 


W ((4 • 8 2 ), q)t-x ,2 = 


A 


(4-8 2 ),3,0,l 


<7-1 


1/3 


(7.3) 


FIG. 4: Plot of the ratios Rt r i,q-,(i,k)/(i,i ) f° r k — 2 to k = 5 as 
functions of q for 4 < q < 6 . From bottom to top, the curves refer 
to k = 2, k = 3, k = 4, and k = 5, respectively. 


where the general expression for D n is given in Eq. (|B1|) . 
Ref. 0 noted that as q increases beyond the lowest values 
above x(^- c ) = 2 , this lower bound rapidly approaches 
the actual value of W(hc,q) (see Table I in 0 ), where 
the latter was determined by a Monte-Carlo simulation 
checked for larger q with a large-q series approximation. 
For example, for example, for q = 3, 4, 5, Rhc,q-, 1,1 is 


where A( 4 .g 2) i3 0 ,i is the largest (real) root of the cubic 
equation (IC4I) in Appendix [C] Correspondingly, W((4- 
8 2 ),y) > W((4- 8 2 ),y)qi, 2 . We analyze the small-y ex¬ 
pansion of W((4 ■ 8 2 ), y)i- 1,2 below. 

VIII. 3 • 6 • 3 • 6 (KAGOME) LATTICE 

In this section we consider the (3 ■ 6 ■ 3 • 6 ) lattice, 
commonly called the kagome lattice (which we shall ab¬ 
breviate as hag). Using the CCM method with 6=1 
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FIG. 5: Section of the (3 • 6 • 3 • 6) (kagome) lattice. 


and k = 1, Ref. {?] derived the lower bound W(kag, q) > 
W(kag,q)i-,x,u where 


7~)2/3 7-)l/3 

W(kag,q) e . AA = ^ _ « 


( 8 . 1 ) 


Equivalently, W(kag,y ) > W(kag, where Q 

W{kag , y ) /;lil = (1 - y 2 ) 2/3 (l + y 5 ) 1 / 3 . (8.2) 


The zigzag path used in the derivation of this lower 
bound was described in detail in Ref. Q • Here, we again 
take 6=1 and k = 1 but use a different type of path. A 
section of the kagome lattice is shown in Fig. [5] Rather 
than the zigzag path used in [7|, we choose the path to be 
given horizontal line in Fig. [5] The matrix T then links 
the proper ^-coloring of the vertices on this line, the ver¬ 
tices between this line and, say, the line above it, and the 
vertices on this higher-lying horizontal line. It turns out 
that the use of this different path yields a slightly more 
restrictive lower bound, which we shall indicate with a 
prime, namely W(kag,q) > W(kag,q)' e . 1 1 , where 


J 


W{kag,qY e . ltl = 


r(<?-2) 


g 4 - 6 q 3 + 14g 2 - 16 q + 10 + 


q 8 - 12<f + 64g 6 - 200g 5 + 404g 4 - 548g 3 + 500g 2 - 292 q + 92 


11/2 


1/3 


(8.3) 


2(g-l) 2 


Equivalently, we have W(kag,y) > W(kag 1 y )' e . 11; where 
W(kag,y)' eil = 


2’ 1/3 (1 + y) 


( l ~y) 


1 — 2y + 2y 2 — 2y 3 + 3y 4 + 


1 — Ay + 8y 2 - 12y 3 + 14 y 4 - 16 y 5 + 16 y 6 - 12 y 7 + 9 y' 


1 / 2 ' 


1/3 


(8.4) 


We find that 

W(kag, q)(-i t i > W(kag,q)i-i,i (8.5) 

The fact that the use of a different path can yield a 
more restrictive bound with the same value of 6 and k was 
already shown for the honeycomb lattice in @,i. Thus, 
both Ref. @ and Ref. [H used the CCM method with 
6=1 and k = 1, but Ref. @ obtained a more restrictive 
lower bound for the honeycomb lattice by using a differ¬ 
ent path. The bounds W(kag, q)e-,i,i and W{kag,q)' i . l ± 
both rapidly approach the actual value of W ( kag , q) as 
q increases beyond the chromatic number, y(fca</) = 3. 
Below we shall show how the slight improvement with 


the new bound is manifested in the respective small-?/ 
expansions of W(kag, q)t-i,i and W(kag, q)' tl In pass¬ 
ing, we note that we have also studied generalizations of 
the CCM method for some other Archimedean lattices. 


IX. sq d LATTICE 

So far, we have considered planar lattices. The coloring 
compatibility matrix method and our generalizations of 
it, also apply to a subclass of nonplanar lattices, namely 
the subclass that can be constructed starting from a pla¬ 
nar lattice and adding edges between vertices on the orig- 
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inal planar lattice. An example of this is the sq d lattice. 
As noted above, the sqd lattice is formed from the square 
lattice by adding edges (bonds) connecting the two sets 
of diagonal next-nearest-neighbor vertices in each square. 
Thus, the vertices and edges in each square form a K 4 
graph. (Here, the Kn graph is the graph with N ver¬ 
tices such that each vertex is connected to every other 
vertex by one edge.) Although an individual K 4 graph 
is planar, the sqd lattice is nonplanar. This lattice has 
coordination number A sgd = 8 and chromatic number 
x( s Qd) = 4. Although it is not 4-partite, an analysis of 
the way in which the number of proper 4-colorings of the 
vertices of a section of the sqd lattice grows with its area 
shows that W(sq d , 4) = 1. 

Using the b = 1, k = 1 CCM, Ref. Q derived the lower 
bound W(sqd,q) > W(sq d , q)e-i,i, where 


W(sq d ,q)t, 1,1 = 


Asq,; ,2,0,1 

Ac,0,1 


(g ~ 2)(g — 3) 
< 2-1 


(9.1) 


A. b = 2, k = 1 


q 4 - 12<f + 54 q 2 - 112g + 97 


11/2 


■ (9-2) 


We thus derive the new lower bound W(sq d ,q) > 
W(sq d ,q)e- 2 , 1 , where W(sq d ,q)e- 2 ,1 = A sgd , 3i0 ,i/A S9s , 2 , 0 , 1 , 
i.e., 


W(sq d ,q)e- 2,1 = 


q 2 — 6 q + 11 + 


q 4 — 12q 3 + 54q 2 — 112q + 97 


11/2 


2(9-2) 


(9.3) 


From these explicit analytic results, we find 

Rsc jd,g;(2,l)/(l,l) A 1 . (9-4) 

That is, our new lower bound W{sq dl q)e, 2,1 is larger and 
hence more restrictive than the one obtained in [6[. 

The corresponding lower bounds for the reduced W 
functions are W(sq d ,y) > W(sq d ,y)e-, 1 , 1 , where 


For our first generalization, namely b = 2 and k = 1, 
we need the dominant A for a cyclic strip of the sq d lattice 
of width L y = 3, which is da 


W(sq d ,y)e-i,i = (1 — 2 /)(l — 2y)(l + y ) 3 

= 1^4y 2 ~2y 3 + 3y 4 + 2 y 5 


A 


sqd, 3,0,1 


(g~3) 

2 


6 q + 11 


(9.5) 


and W{sqd, y) > W(sq d ,y)e- 2 . 1 , where 


(i + yf 


W(sqd,y)t- 2,1 = 


1 - 4y + 6 y 2 + 


1-8 y + 24 y 2 - 36 y 3 + 28 y 4 


1 1/2 


2(1 — y) 


(9.6) 


B. b = 1, k = 2 


For b = l and k = 2, we derive the lower bound W(sqd,q) > W(sqd,q)e- 1 , 2 , where 


W(sq d ,q)t, 1,2 


As 9d ,3,0,l 

1/2 

'(g- 3) 

q 2 - 6 q + 11 + 

V - 12 q 3 + 54 9 2 - 112g + 97^ 

!/2- _ 

- Ac, 0,1 ■ 


2 (g — 1) 



(9.7) 


Equivalently, W{sq d ,y ) > W{sq d , y)t, 1 , 2 , where 


W{sq d ,y)t- 1,2 = -j={l + y ) 3 


(1 - 2y) [l - Ay + 6 y 2 + (l - 8 y + 24 y 2 - 36 y 3 + 28 y 4 ) 


1/21 


1/2 


(9.8) 


mentary to, that of a Taylor series expansion, in this 
case, a small-y expansion. The lower bound is valid for 
any value of q that is physical, but need not, a priori , 
be an accurate approximation to the actual function. In 

A lower bound on a function such as W(A, q) or 
W (A, y) plays a role that is different from, and comple- 


X. SMALL-y EXPANSIONS OF NEW LOWER 
BOUNDS 

A. General 































13 


contrast, the large-g (equivalently, small-y) Taylor series 
expansion is an approximation to the function itself and, 
within its radius of convergence, it satisfies the usual Tay¬ 
lor series convergence properties. Thus, if one truncates 
this series to a fixed order of expansion, then it becomes 
a progressively more accurate approximate as the expan¬ 
sion variable becomes smaller, and for a fixed value of 
the expansion variable, it becomes a more accurate ex¬ 
pansion as one includes more terms. 

A lower bound on a function W(A, y) need not, a pri¬ 
ori , agree with the terms in the small-y Taylor series 
expansion of this function. Some explicit examples of 
this are given in Appendix [A] Interestingly, as discussed 
in 0-0, the lower bounds derived there do agree with 
these small-y series to a number of orders in y (listed for 
Archimedean lattices in Table III and for the duals of 
Archimedean lattices in Table IV of Ref. 0 )- 

It is thus clearly of interest to carry out a similar com¬ 
parison to determine the extent to which our new lower 
bounds, which we have shown improve upon those in [§} 
and 1-0, agree with the respective small-y expansions 
to higher order. We do this in the present section, show¬ 
ing that our new lower bounds are not only more strin¬ 
gent than the earlier ones, but also agree with the small- 
y expansions of W (A, y) to higher order in y than these 
earlier lower bounds. 

Because W(A, y)t-b,k is a lower bound on W( A, y), one 
can draw one immediate inference concerning the com¬ 
parison of the small-y Taylor series for these two func¬ 
tions, namely that for a given lattice A, if the small-y 
Taylor series of W(A, y)e-b,k coincides with the small-y 
series for W(A,y) to order 0(y lc ), inclusive, then the 
difference 

W(A,y) - W (A, y)e-,b,k = m c +iy lc+1 , with m c+ i > 0 . 

_ ( 10 . 1 ) 

Thus, for example, with the 0(y v+1 ) term in W( A, y) de¬ 
noted t0A,i o +i and with the 0(y v+1 ) term in W(A, y)e-,b,k 
denoted u>A;6,fc;i c +i, we have 

WA;i„+l > W A ;b,fc; ic + 1 • (10.2) 

We discuss a subtlety in this comparison. One should 
first show that the small-y expansion is, in fact, a Taylor 
series expansion, i.e., that W(A,y) is an analytic func¬ 
tion at y = 0 in the complex y plane, or equivalently, that 
W r (A,q) is an analytic function at 1/q = 0 in the com¬ 
plex plane of the variable 1/q. In fact, there are families 
of N- vertex graphs Gn such that W r ({G}, q) is not ana¬ 
lytic at 1/q = 0 [28], where here {G} denotes the formal 
limit limAr-»oo Gn- This is a consequence of the property 
that the accumulation set of zeros of the chromatic poly¬ 
nomial P{Gn , q), denoted B , extends to infinite |g| in the 
q plane, or equivalently, to the point 1/q = 0 in the 1/q 
plane. (The zeros of P(G, q) are denoted as the chromatic 
zeros of G.) Refs. [28[ constructed and analyzed various 
families of graphs for which this is the case. For regular 
(vertex-transitive) A-vertex graphs Ga,at of a lattice A 
with either free or periodic (or twisted periodic) bound¬ 
ary conditions, the resultant W r (A, q) functions obtained 


in the N oo limit are analytic at 1/q = 0. This fol¬ 
lows because a necessary condition that B extends to 
infinitely large |g| as N — > oo is that the chromatic zeros 
of Ga,at have magnitudes |g| —» oo in this limit. How¬ 
ever, a vertex-transitive graph G has the property that 
all vertices have the same degree, A and a chromatic 
zero of G has a magnitude bounded above as |g| < 8.4A 
[29l ]. So for the N — > 00 limit of a regular lattice graph A, 
W r ( A, q) is analytic at 1/q = 0 and equivalently, W(A, y) 
is analytic at y = 0 , and the corresponding series expan¬ 
sions in powers of 1/q and powers of y are Taylor series 
expansions. 

B. Square Lattice 

The small-y expansion of W ( sq , y) is [lb] 

W(s qi y) = 1 + y 3 + y 7 + 3y 8 + 4y 9 + 3y 10 

+ 3y n + 0(y 12 ) . (10.3) 

This series and several others for regular lattices are 
known to higher order than we list; we only display the 
various series up to the respective orders that are rel¬ 
evant for the comparison with our lower bounds. As 
is evident from Eq. (14.31) . the previous lower bound 
W(sq,y) e - 1,1 = l + y 3 a coincides with the small-y series 
to 0 (y 6 ), inclusive. 

We list below the small-y expansions of the various new 
lower bound functions W(sq, y)t,b,k that we have derived 
with b > 2 and k = 1 and with b = 1 , k > 2 : 

W ( sq , y) e - 2,1 = 1 + y 3 + y 7 + 3y 8 + 3y 9 + 0(y 10 ) 

(10.4) 

W(sq,y)e- 1,2 = l + y 3 + ^y 7 + | y s + ^ y 9 + 0 (y 10 ) 

(10.5) 

and 

W ( sq , y) e . 1,3 = 1 + y 3 + ^y 7 + 2y 8 + ^y 9 + 0(y 10 ) . 

( 10 . 6 ) 

Comparing the small-y expansion of our new lower 
bound function W(sq,y)e-b, 1 with b = 2 , as well as 
the old lower bound function W(sq, y)qi,i, with the ac¬ 
tual small-y series for W(sq,y) in Eq. (110.31) . we can 
make several observations. First, the small-y expansions 
for W(sq, y)t, 2 ,\ coincides with the small-y expansion of 
W(sq, y) to 0(y s ), inclusive, which is an improvement by 
two orders in powers of y as compared with W(sq, y)t, 1,1 
(see Eq. m)- Since increasing b (with k fixed) im¬ 
proves the accuracy of the lower bound, it follows that 
W(sq, y)e-b,i will also coincide with the series for W(sq, y) 
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to at least 0(y 8 ) for 6 > 3 as well as for 6 = 2. Moreover, 
although the respective coefficients of y 9 in the series for 
W(sq,y)i-i,i and W(sq, y)t, 2 ,i, namely 0 and 3, do not 
match the coefficient of y 9 in the actual small-?/ expansion 
of W(sq,y), which is 4, one can see that as b increases 
from 1 to 2, this coefficient of the y 9 term increases to¬ 
ward the exact coefficient. 

Regarding the matching of terms in the small-?/ expan¬ 
sions of the W(sq , y)t-b, i, as compared with W ( sq , y)e-i,ki 
that we have calculated, we find that this matching is bet¬ 
ter by two orders for the W(sq, y)e-,b, l than W(sq, y)t,i,k- 
That is, for the k values that we have calculated, namely 
k = 2, 3, the lower bounds W(sq, y)t,i,k match the small- 
y expansion of W(sq,y ) to order 0(y 6 ), the same order 
as W(sq,y)i-iA- 

A related property of our lower bounds for a general 
lattice A and, in particular, for the square lattice, follows 
as a consequence of the theorem (13.91) and (13.101) : with 
6=1, since the lower bound W(A, y)e-,i,k is a monoton- 
ically increasing function of k, the degree of matching 
of coefficients in the small-?/ expansion for W(A, y) must 
improve monotonically as k is increased. A priori, this 
improvement could be manifested in two ways (or a com¬ 
bination of the two): (i) as k is increased, coefficients of 
terms of higher order in y are exactly matched, or (ii) 
the coefficient of a given term of a certain order in y ap¬ 
proaches monotonically toward the exact value. For the 
present lattice A = sq, we see that, for the W(sq,q)i-i^ 
that we have calculated, the latter type of behavior, (ii), 
occurs. That is, as we increase k from 1 to 2 to 3, the coef¬ 
ficient of the y‘ term in the small-?/ series for W(sq, y)e-,i,k 
increases from 0 to 1/2 to 2/3, moving toward the exact 
value of 1. This is similar to the behavior that we ob¬ 
served with the respective coefficients of the y 9 term in 
the small-?/ expansions of W(sq, y)t,b,i as compared with 
the exact value. This type of behavior is in accord with 
the inequality (110.21) . 

Regarding the relative ordering of the various lower 
bounds that we have obtained, from the small-?/ expan¬ 
sion, we find, for large q, the ordering 

W(sq,y) > W(sq,y)e- 3,1 > W(sq,y)c- 2,1 
> W(sq,y)e- 1 , 3 > W(sq,y)i- 1 , 2 > W(sq,y)t-i,i ■ 

(10.7) 

In fact, we find that this ordering also extends down 
to the lowest value where we apply our lower bounds, 
namely q = 3. For bounds on W(sq, 4) and W(sq, 5), see 

& ' 

C. Triangular Lattice 

The small-?/ expansion of W(tri,y) is 0 

W(tri, y) = 1 — 2 y 2 + y 4 + y 5 + 5 y 6 + 16 y 7 + 47 y 8 


+ 134 y 9 + 0(y 10 ) . (10.8) 

As is evident from Eq. (15.21) . the previous lower bound 
W(tri,y)e-i,i = (1 — y 2 ) 2 fa, [3/ matches the small-?/ series 
to 0(y 4 ), inclusive. 

We list below the small-?/ expansions of the various new 
lower bounds W(tri, y)e-b,k that we have derived with 
6 > 2 and k = 1, and with 6 = 1, k > 2: 

W(tri, y)e-, 2 ,i =1-2 y 2 + y 4 + y 5 + 5 y 6 + 14 y 7 + 0(y s ) 

(10.9) 

W(tri, ?/)/? ; i ,2 = 1 - 2y 2 + y 4 + i y 5 + ^ y 6 + 0(y 7 ) 

( 10 . 10 ) 

and 

_ 9 10 

w(tri, y)t; 1,3 = 1 - 2?/ 2 + y 4 + -y 5 + —y 9 + 0(y 7 ) . 

( 10 . 11 ) 

Comparing these with the small-?/ series for W(tri, ?/)qi,i, 
we find that, among (110.91) - (IIP. 1 II) . the greatest match¬ 
ing of terms is achieved with (110.91) . i.e., by increasing 
6. Specifically, the small-?/ expansion for W(tri, y)e- 2 ,i 
matches the small-?/ expansion of W(tri,y ) to 0(y 6 ) in¬ 
clusive, which is an improvement by two orders in y as 
compared with W(tri, y)e-,i,i- This increase by two orders 
in y is the same amount of improvement that we found 
for our lower bound for the square lattice, W(sq,y)t, 2,1 
as compared with W(sq, y)t-i,i. 

As was true of the lower bounds for the square lat¬ 
tice, the lower bounds W(tri, y)t,i,k with k = 2 and 
k = 3 coincide with the small-?/ series for W(tri,y) to 
the same order, namely 0(y 4 ), as W(tri, y)qi,i. How¬ 
ever, as k increases from 1 to 2 to 3, the coefficient of 
the first unmatched term in the respective small-?/ se¬ 
ries for W(tri,y)e- t i,k, viz., the y 5 term, increases from 0 
to 1/2 to 2/3, moving toward the exact value of 1. An 
inequality that follows from the theorem (13.91) and gen¬ 
eral result (13.101) . is that with 6=1, W(tri,y)t,\,k is a 
monotonically increasing function of k. 

Concerning the relative ordering of the various lower 
bounds that we have obtained, from the small-?/ expan¬ 
sion, we find, for large q, the ordering 

W(tri,y) > W(tri,y)e- 2 ,i > W(tri,y)e- 1,3 > 

> W(tri,y)e-,i ,2 > W(tri,y)e ; i,i ■ (10.12) 

Indeed, we find that this ordering also extends down to 
the lowest value where we apply our bounds, namely q = 

4. 
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D. Honeycomb Lattice 


The small-y expansion of W (he, y) is [2C| 
W(hc,y ) = 1 + + »" + 2 # 

+ y u - + 0(/«). 

The previous lower bound W(hc, y)i-i,i = 
Q-0] has the small -y expansion 


+ -y 
2 y 

(10.13) 
(1 + y 5 ) 1 / 2 


W(he, y)e-i,i = 1 + \y 3 - ± y 10 + ±y 15 + 0 (y 20 ) . 

(10.14) 

Thus, as was noted in 1-0 , this small-?/ expansion coin¬ 
cides with the small-?/ expansion of W(he, y) to the quite 
high order 0 (y w ). 

We list below the small-?/ expansions of the various new 
lower bound functions W(he, y)e-b,k that we have derived 
with b > 2 and k = 1 and with b = 1, k > 2: 


W(hc,y) e . 2,1 = l + i?/ 5 -l?/ 10 + ?/ 11 + 2?/ 12 + ?/ 13 + O(?/ 15 ) 

(10.15) 

and 

W(he, y)i- 1,2 = 1 + ^y 5 - ^y 10 + ^y 11 + y 12 + 0(y 13 ) . 

(10.16) 

As with the square and triangular lattices, we find that 
among (U0.15I) - (I10.16I) . the greatest matching of terms 
is achieved with (110.151) . i.e., by increasing b. Specifi¬ 
cally, the small-?/ expansion for W(hc, y)e ; 2,1 matches the 
small-?/ expansion of W(hc,y) to 0(y 12 ) inclusive, which 
is an improvement by two orders in y as compared with 
W(hc, y)e-i,i. 

The theorem (l3fil) and corollary (13.101) imply that 
W(hc, > W(hc, q)t,i,i, and this inequality is re¬ 

flected in the degree of matching of the small-?/ expan¬ 
sions for the corresponding functions W(hc,y)t,i ,2 and 
W(hc,y)t,i,i- Although W(hc, y)t-, 1,2 does not increase 
the order of matching, as compared with W(hc, y)i-i t i, 
it begins the process of building up a nonzero coeffi¬ 
cient for a y 11 term, which was zero in the expansion 
of W(hc, y)e-i,i- Specifically, the small-?/ expansion of 
W(hc,y)t, 1,2 contains a y 11 term with coefficient 1/2, 
building toward the exact coefficient, 1, of y 11 in (110.131) . 


E. 4 • 8 J Lattice 

We next consider a (bipartite) heteropolygonal 
Archimedean lattice, namely the (4 • 8 2 ) lattice. The 
small-?/ expansion of W(( 4 • 8 2 ), y) is !6f, 0 

W((4-8 2 ), y ) = 1 + \y 3 - ^y 6 + \ y 7 + ^y 9 


+ ^y 10 -^y 12 + O(y 13 ) ■ (10.17) 
The small-?/ expansion of the lower bound obtained in 

@,0,W((4-8 2 ), ?/),;!,!, is 

w ((4 ■ 8 2 ), y)t;i,i = 1 + \y 3 - -^y 6 + \v 7 + ^y 9 

1 10 77 12 3 13 3 14 231 15 16\ 

+ y -yiv -yv -55 y + ^ y +°(y )■ 

(10.18) 

As was noted in @,0 , this coincides with the small-?/ 
expansion of W((4-8 2 ), y) to the quite high order 0(y 12 ). 

We list below the small-?/ expansions of the various 
new lower bound functions W((4-8 2 ), y)e-,b,k that we have 
derived with b > 2 and k = 1 and with b = 1, k > 2: 

W((4-8 2 ),y) e;2 ,i = l + ^-y 3 - J-y 6 + jy 7 

7 9 1 10 77 12 189 13 14 . 

+ 27 y +34 y - 211 y + ~ w v + ) 

(10.19) 


and 


W(( 4 • 8 2 ), y)t ,\,2 = 1 + \y 3 - ^y 6 + \y 7 

7 9 1 10 77 12 9 3 13 45 14 15 

+ 2?y+^y -guy + ^y + 32 y + °( y )• 

( 10 . 20 ) 

Evidently, the small-y series expansions of W((4 • 
8 2 ),y)t; 2 ,i and W((4 • 8 2 ),y)t,i ,2 match the small-y ex¬ 
pansion of W((4 • 8 2 ),y) to at least the same order as 
W((4 ■ 8 2 ),y)e ; i i i. Further, we observe that for small-y, 

W((4 • 8 2 ),y)e- 2 ,i > W((4 • 8 2 ), y) t; 1>2 > W((4 ■ 8 2 ), y )* ;1)1 

( 10 . 21 ) 

F. 3 ■ 6 • 3 ■ 6 (kagome) Lattice 
The small-y expansion of W(kag,y) is [Tf 

W(kag, y) = 1 - |y 2 - ^y 4 + iy 5 

-^y 6 -^y 7 -^y s + 0(y 9 ). ( 10 . 22 ) 

As was discussed in [3] , the small-y expansion of the b = 
1 , k = 1 lower bound W(kag,y)e ; 1,1 derived there (listed 
above as Eq. (18.11) 1 coincides to 0(y 8 ) with the small-y 
series for the actual quantity W(kag,y). Explicitly, 

W(kag, y)e-i,i = 1 - ^y 2 ~ ^y 4 + ^ 2/ 5 
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4 6 2 7 7 


1 9 95 1 0 


34 y 32 y 35 y 33 y 36 y 
+0(y n ) . (10.23) 

Our new bound has the small-y expansion 

W ( kag , y)' e . 1}1 = 1 -^y 2 - ^y 4 + ^y 5 


4 ..6 2 - 7 


® „9 , ^,10 


34 y 32 f 35 y 33 ; 

+0{y n ) . 


3 6 


(10.24) 


Thus, 

W{kag, y)' e . 1A - W{kag, y)* ; i,i = y y 9 + 0(y 10 ) • 

(10.25) 

One could derive similar lower bounds for other 
Archimedean lattices not considered here, e.g., the 3 • 
12-12 lattice [1U. 


G. sqd Lattice 


deriving lower bounds on the ground-state degeneracy 
per site, W{A,q), of the y-state Potts antiferromagnet 
on several different lattices A. Our first generalization is 
to consider a coloring compatibility matrix that relates a 
strip of width 6 > 2 vertices to an adjacent strip of the 
same width. Our second generalization is to consider a 
coloring compatibility matrix that acts k > 2 times in re¬ 
lating a path on A to an adjacent parallel path. We have 
applied these generalizations to obtain new lower bounds 
on W{A,q ), denoted W{A,q)t,b,k- In this notation, the 
lower bounds previously derived in i-i have 6=1 and 
k = 1. One of the interesting properties of these bounds 
W (A, q)t,x, 1 obtained in was that as q increases be¬ 

yond x(A) they rapidly approach quite close to the actual 
respective values of W{A,q). We have shown that our 
new lower bounds are slightly more restrictive than these 
previous lower bounds, and consequently are even closer 
to the actual values W (A, q). We have demonstrated how 
this is manifested in the matching to higher-order terms 
with the large-y (small-y) Taylor series expansions for the 
corresponding functions W{A,y) for the various lattices 
that we have considered. 


Since the lower bound W{sq d ,y)t-,i,i derived in [|| and 
given above in Eq. (19.511 is a polynomial, it is identical 
to its small-y Taylor series expansion. 

Expanding W{sq d ,y)e- 2 , 1 , we find 

w{sq d , y)i; 2 ,i = 1-4y 2 - y 3 + 6y 4 + 5y 5 + 0(y 6 ) . 
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(10.26) 

Similarly, 

W(sq d ,y)e-,1,2 = 1- 4y 2 - ^y 3 + |y 4 + ^y 5 + 0 (y 6 ) 

(10.27) 

bE (sq d , y) e - 1,3 = 1 - 4y 2 - ^y 3 + 5y 4 + y y 5 + 0(y 6 ) . 

(10.28) 

From these expansions we find, for large q, the ordering 

W{sqd, y) > W{sq d ,y)i- 2 ,1 > W{sq d ,y)t- 1,2 > W{sq d ,y)i- 1,1 

(10.29) 

This is the same ordering that we found for the other 
lattices. 

XI. CONCLUSIONS 

Nonzero ground-state entropy per site, So, and the as¬ 
sociated ground-state degeneracy per site, W = e s °^ kB , 
are of fundamental importance in statistical mechanics. 

In this paper we have presented generalized methods for 


Appendix A: W{A rp .,q) at q = x(A) 

We mention here a subtlety that results from the non¬ 
commutativity in the limits (12.121) . An r-partite rp. 
graph with N vertices, G rp , t N has chromatic number 
x(Crp. I Jv) = t. One equivalent definition of an r-partite 
graph is that its chromatic polynomial, evaluated at 
q = r, satisfies 

P{G rp .^ N ,r) = r\ (Al) 

The square and honeycomb lattices are bipartite (as are 
the (4-6-12), and (4-8-8) lattices, among Archimedean 
lattices), while the triangular lattice is tripartite (for oth¬ 
ers Archimedean lattices and their planar duals, see, e.g., 
Tables I and II in 0)- It follows that, with the D^q defi¬ 
nition for W{ A, q), namely setting q = r and then taking 
the N — > 00 limit in Eq. & one has 

W(A rp .,r) = l. (A2) 

As discussed in 0, because of the noncommutativity 
(12.121) . if instead of setting q = r, evaluating P{G rp .,r), 
and then taking the N —> 00 , one first takes N —> 00 with 
q in the vicinity of r, and then performs the limit q —> r, 
one can, in general, get a different result for W{A,q). 
Indeed, this is the case for many lattice strips of regular 
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lattices of a fixed width L y , an arbitrary length, L x , and 
various transverse and longitudinal boundary conditions 
BEJEJil. The coloring problem on a given lattice 
A is of interest for q > x(A), since this is the minimum 
(integer) value of q for which one can carry out a proper 
(/-coloring of the vertices of A. In a number of cases, 
x(A) < q c { A). If one considers W(A,q) for q < q c { A), 
then one must deal with the generic noncommutativity 
in the limits (12.1211 fl3 ). Here we always use the order 
D^q, i.e., we fix q to a given value and then take N —> oo. 
Actually, in view of the results m and (1A2I) . for the 
square and honeycomb lattices, W(sq,2) = W(hc,2) = 
1, and for the triangular lattice, W{tri, 3) = 1. Since 
our new lower bounds are intended for practical use and 
since one already knows (with the Djvq definition) the 
values of W(sq, 2), W{hc,2), and W{tri,3) exactly, we 
may restrict our analysis to the application of our new 
bounds in the range q > 3 for the square and honeycomb 
lattices and to the range q > 4 for the triangular lattice. 

For reference, we recall an elementary lower bound on 
P{G, q) and hence on liniAr_ ) . 00 P{G, q) r ^ N , where G is an 
N-vertex graph. If G is bipartite {bp.) then one can assign 
a color to all of the vertices of the even subgraph in any 
of q ways and then one can assign one of the remaining 
q — 1 colors to each of the vertices on the odd subgraph 
independently, so P{Gb P .,q ) > q(q — 1 ) N ^ 2 . Hence, for 
a bipartite lattice, denoting Ab p . as the N ^ oo limit 
of Gbp., one has W{Ab p .,q) > (q — l) 1 / 2 . Both of these 
lower bounds are realized as equalities only in the case 
q = 2. More generally, if G rp , is an r-partite graph and 
A xp . — limjv —yoo Arp. ,jv •> then 


r -2 


P(G rp .,q)> [n(9- a )][9-( r -l)]* /r ( A3 ) 

s—0 

and hence 

W(A rp .,q)> [q — {r — l)] 1 /' r . (A4) 


Thus, for example, one has the elementary lower bounds 
W{sq, q) > {q — l) 1 / 2 and W ( tri, q) > {q — 2) 1 / 3 , etc. 

For q > r on A rp ., the lower bound (IA4D is less strin¬ 
gent than the ones derived in i-i and here via coloring 
matrix methods. Indeed, these lower bounds illustrate 
the fact noted in the text, namely that, a priori , a lower 
bound need not agree with terms in the large-g expan¬ 
sion of W r {A,q) or the equivalent small-// expansion of 
W (A, y ). For example, for the square and honeycomb 
lattices, the r = 2 special cases of (IA4I) read, for q > 2, 


W{sq,q)>{q- 1) 1/2 

(A5) 

W{hc,q)> {q-l) 1/2 . 

(A6) 


value) as q increases above 2. The corresponding lower 
bounds in terms of W{sq, y) and W{he, y ) are 


W{sq,y) > (1 + y)\fy 

(A7) 

and 


W{hc, y) > y/y(l + y) . 

(AS) 

Rather than matching any terms in the respective small- 
v exDansions (110.31) and (110.131). the right-hand sides of 
these lower bounds vanish for small y. Similarly, since 

the triangular lattice is tripartite, the r = 
of (IA4I) yields the lower bound, for q > 3, 

3 special case 

W{tri,q) > {q- 2) 1/3 . 

(A9) 

In terms of W{tri,y), this is 


W{tri,y) > y 2/3 {l - y) 1/3 { 1 + y) 2 

(A10) 


Again, for small y, this vanishes rather than matching 
any of the terms of the small-// expansion (110.81) . Thus, as 
noted, a lower bound need not match any of the terms in 
the small-// expansion. This emphasizes how impressive 
the new lower bounds are in their matching of these terms 
in the small-// expansions for the various lattices to high 
order. 


Appendix B: Lower Bounds W{A, and 

W (A, //)(;!, i for Archimedean Lattices 


A number of general results were proved in Ref. 
concerning lower bounds which, in the notation of this 
paper, are W{A, q)e ; i,i and W{A, y)e-,i,i- These results 
in [zj applied for all eleven Archimedean lattices. We 
have given the notation for an Archimedean lattice in 
the text. These lower bounds from |7] are relevant here 
because we compare our new lower bounds W (A, q)e-b,k 
and the corresponding lower bounds W{A, y)(.-,b,k with 
b > 2 and/or k > 2 to these earlier ones with b = k = 1. 
(Ref. [7] also gave lower bounds for the planar duals of 
the Archimedean lattices; we do not list these here but 
instead refer the reader to [7].) 

The chromatic polynomial of a circuit graph is 
P{C n ,q) = {q — l) n + {q — 1)(—l) n . Since this chro¬ 
matic polynomial has q{q — 1) as a factor, we can write 
it as P{C n ,q) = q{q - 1 )D n {q), where 


D n {q) 


P{C n ,q) 

q{q~l) 




s =0 


n — 1 
s 


jn—2—s 


(Bl) 


Ref. Q proved the following general lower bounds for 
an Archimedean lattice, A = 01 ■ p“ 4 ) (where we add the 
subscripts 1,1 to indicate 6=1 and k = 1 to match our 
current notation for W{A,q)^b,k)'- 


Since W{A,q) ~ q for large q, these lower bounds be- W (TTp® 4 ), q) > W( (TTp“‘)> <7j , (B2) 

comes progressively worse (i.e., farther from the actual ^ ^ i ' ^L 1 
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where 


W 


Ikrul 


q~ 1 


A summary of these for Archimedean lattices is given in 
Table IV of 0. 

(B3) 


Here, the {*} in the product label the set of p^-gons in¬ 
volved in A and v p (with p = pi here) was defined in Eq. 
(2.10) of 0 . 

This lower bound takes a somewhat simpler form in 
terms of the related function W(A,y)g, namely, 

w(([[pr),y) > , (B4) 


Appendix C: Higher-Degree Algebraic Equations for 
Certain Aa,z,„,o,i 


where 


W 


Ibr)>v\ 


n[i + (-i ) Pi y Pi -\ 


(B5) 


In this appendix we list some algebraic equations of 
degree higher than 2 that are used in the text. The cubic 
equation whose largest (real) root is A sq , 4 , 0 , 1 , used for our 
lower bound W(sq, q)t, 1 , 3 , is 


A 3 - ( q 4 - 7 q 3 + 23 q 2 - 41 q + 33) A 1 2 + (2 q 6 - 23 q 5 + 116g 4 - 329 q 3 + 553 q 2 - 517 q + 207) A 

- q 8 + 16 q 7 - 112g 6 + 449g 5 - 1130g 4 + 1829g 3 - 1858g 2 + 1084g - 279 = 0 . (Cl) 

The quartic equation whose largest (real) root is Xtri, 4 , 0 , 1 , used for our lower bound W(tri, q)e-i,a, is 
A 4 - (g 4 - 10 q 3 + 42 q 2 - 88 q + 76) A 3 + (g - 2)(g - 3) 2 (3 q 3 - 22 q 2 + 60 q - 60) A 2 
- (q — 2 ) 2 (q — 3) 3 (3q 3 — 21 g 2 + 51g — 43) \+ (q — 2 ) 6 (q — 3 ) 4 = 0 . (C2) 

The cubic equation whose largest (real) root is A^c^^p used in our lower bound W(hc, < 7 ) 42,1 is 

A 3 - (q 6 - 8 q 5 + 28 q 4 - 56 q 3 + 71 q 2 - 58 q + 26) A 2 
+ (q - l) 2 (q 6 - 10 q 5 + 43 q 4 - 102 q 3 + 144g 2 - 120 g + 49) A - (q - 1 ) 4 (q - 2 ) 2 = 0 . (C3) 

The cubic equation whose largest (real) root is A( 4 . 82 ), 3 ,o,i, used in our bound W((4 • 8 2 ),g)i ; i i 2 is 

A 3 - (g 12 - 16g n + 120g 10 - 558 q 9 + 1794 q 8 - 4212g 7 + 7437g 6 - 10018g 5 + 10324g 4 
— 8064g 3 + 4648g 2 - 1854g + 414) A 2 

+ {q - l) 4 (q 12 - 20g 11 + 188g 10 - 1094g 9 + 4375g 8 - 12640g 7 + 27033g 6 - 43164g 5 
+51235g 4 - 44380g 3 + 26931g 2 - 10462g + 2017) A 

-(g - 1) 8 (9 - 2) 2 (g - 3) 2 (g 2 - 3g + 3) 2 = 0 . (C4) 
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TABLE I: Values of R sq ,q;e;b,k f° r 9 = 3 and some illustrative 
values of b and k. 


b 

k 

W(sq,3)e- b}k 


1 

1 

1.500000 

0.974279 

2 

1 

1.520518 

0.987605 

3 

1 

1.530340 

0.993985 

1 

2 

1.510224 

0.980919 

1 

3 

1.5162645 

0.984843 

1 

4 

1.520249 

0.987430 

1 

5 

1.523073 

0.989265 


TABLE II: Values of Rtri,q;t;b,k f° r 9 = 1 and some illustrative 
values of b and k. 


b 

k 

W{tri,A)i- b ,k 

Rtri,4-,£-,b,k 

1 

1 

1.333333 

0.912618 

2 

1 

1.390388 

0.951670 

3 

1 

1.427052 

0.976765 

1 

2 

1.361562 

0.931939 

1 

3 

1.380569 

0.944949 

1 

4 

1.393923 

0.954089 

1 

5 

1.403672 

0.960762 




